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Efficient  Solution  For  Electromagnetic  Scattering 
Using  the  Dual-Surface  Magnetic-Field  Integral 
Equation  For  Bodies  of  Revolution 


1  INTRODUCTION 

Conventional  solutions  of  the  magnetic-field  integral  equation  (MFIE)  fail  to  produce 
a  unique  solution  at  frequencies  equal  to  the  resonant  frequencies  of  the  interior  cavity.^’^ 
These  spurious  resonances  severely  corrupt  the  numerical  solution  of  the  MFIE.  Yaghjian^ 
proved  that  the  original  MFIE  allows  spurious  resonances  at  these  frequencies  because 
it  does  not  restrict  the  tangential  electric  field  to  zero  on  the  surface  of  the  scatterer. 
Tobin,  Yaghjian,  and  Bell(Woodworth)^  derived  the  dual-surface  magnetic-field  integral 
equation  (DMFIE)  and  applied  it  to  a  three-dimensional,  multi- wavelength,  perfectly 
conducting  body.  This  numerical  solution  showed  that  the  DMFIE  did  indeed  eliminate  the 
spurious  resonances  associated  with  the  interior  cavity  modes.  The  proof  of  uniqueness  and 
derivation  for  both  the  dual-surface  electric-field  integral  equation  (DEFIE)  and  magnetic- 
field  integral  equation  appeared  in  Woodworth  and  Yaghjian.® 

Mautz  and  Harrington®  used  the  combined-field  integral  equation  solution  (CFIE)  to 
eliminate  the  spurious  resonances  associated  with  the  interior  cavity  modes.  The  CFIE, 
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shaped,  perfectly  conducting  bodies,  Digest  of  the  National  Radio  Science  Meeting,  (URSI),  Boulder  CO 

®  Woodworth,  M.B.,  and  Yaghjian,  A.D.  (1991)  Derivation,  application  and  conjugate  gradient  solution  of  the 
dual-surface  integral  equations  for  three-dimensional,  multi- wavelength  perfect  conductors,  in  PIERS- 5:  Applications  of 
the  Conjugate  Gradient  Method  to  Electromagnetic  and  Signal  Analysis,  Sarkar,  T.K.  and  Kong,  J.A.,  ^s.,  Elsevier, 
ch.  4.  Also,  J.  Opt.  Soc.  Am.  A,  11,  April  1994. 

®  Mautz,  J.R.  and  Harrington,  R.F.  (1987)  H-field,  E-field,  and  combined-field  solutions  for  bodies  of  revolution. 
Arch.  Elecktron.  Ubertragungstech.  (Electron.  Commun.),  32(4):157-164. 
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as  the  name  implies,  combines  the  solution  of  the  MFIE  with  the  electric-field  integral 
equation  (EFIE).  Because  the  CFIE  involves  the  MFIE  and  the  EFIE  operators,  additional 
programing  ability  and  computer  run-time  are  required  than  with  the  original  MFIE. 

One  application  of  the  CFIE  and  DMFIE  solutions  is  plane  wave  scattering  from  a 
body  of  revolution.  As  an  alternative  to  the  CFIE,  a  body  of  revolution  solution  was 
formulated  using  the  DMFIE  to  determine  if  it  was  advantageous  to  use  because  of 
increased  accuracy  or  decreased  computer  run-time.^  The  body  of  revolution  formulation  of 
the  DMFIE  solution  parallels  the  body  of  revolution  formulation  in  Mautz  and  Harrington. 
However,  the  numerical  solution  to  the  body  of  revolution  DMFIE  uses  a  single  pulse 
and  impulse  approximation  for  basis  and  testing  functions  instead  of  the  triangular  basis 
and  testing  functions  used  in  Mautz  and  Harrington.®  The  CFIE  solution  requires  a 
more  complicated  set  of  basis  and  testing  functions  because  it  includes  the  EFIE.  The 
different  basis  and  testing  functions  are  necessary  to  accurately  model  the  derivative  of 
the  surface  current  in  the  EFIE  which  is  not  present  in  the  MFIE.  Also,  summations 
for  the  integration  use  a  simple  rectangular  rule  rather  than  the  more  complex  Gaussian 
quadrature  integration  scheme. 

The  purpose  of  this  report  is  to  present  the  development  of  the  DMFIE  for  a  perfectly 
conducting  body  of  revolution  and  its  capability  for  solving  electrically  large  bodies  using 
the  fast  Fourier  transform  (FFT)  and  conjugate  gradient  (CG)  methods.  Also  presented 
is  the  implementation  of  the  dominant  tip  current  variation  required  to  solve  bodies 
with  sharp  tip-like  structures  (such  as  conespheres  with  half-cone  angles  less  than  20°) 
using  the  MFIE  and  DMFIE.  Results  are  presented  with  and  without  the  FFT  and 
CG  methods.  Comparisons  of  monostatic  radar  cross  section  for  the  original  DMFIE, 
the  DMFIE  with  modified  dominant  tip  current  implementation,  and  measurements  are 
presented  for  various  scatterers  with  tips. 


2  STATEMENT  OF  PROBLEM 

It  is  desired  to  determine  the  scattered  far  field  from  a  perfectly  conducting  body  of 
revolution  excited  by  a  plane  wave.  Figure  1  shows  the  geometry  and  the  coordinate  system 
used  for  the  body  of  revolution,  where  p,  0,  and  z  are  the  cylindrical  coordinates  and  t 
and  <j)  form  an  orthogonal  curvilinear  coordinate  system  on  the  body  of  revolution  surface 
S.  Ut  and  are  orthogonal  unit  vectors  in  the  t  and  (j)  directions,  respectively.  Figure 
2  shows  the  coordinate  system  for  plane  wave  scattering,  h  is  the  propagation  vector  for 
the  incident  plane  wave,  $t  defines  the  transmit  angle,  9r  ai^d  4>r  are  the  receiver  spherical 
angles  at  which  the  scattered  far  field  is  observed,  and  h  ^  the  receiver  propagation 
vector,  which  points  from  the  receiver  location  to  the  origin,  kt  is  limited  to  the  xz-plane 
=:  0).  Figure  2  shows  the  unit  vectors  ttg,  and  in  the  Of,  y,  Or,  and  (f)r  directions 

respectively.  The  notation  follows  basically  that  of  Mautz  and  Harrington. 

Considered  separately  is  an  incident  plane  wave  (H^)  defined  by 

^hmitz,  J.L.  (1991)Dual-Surface  Magnetic- Field  Integral  Solution  for  Bodies  of  Revolution,  RLrTR-91-139, 
Hanscom  AFB,  MA,  ADA260725. 
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Figure  1.  Body  of  Revolution  and  Coordinate  System. 


=  (1) 

and  an  incident  plane  wave  defined  by 

Hinc  =  {-ktXul)e-^^^-^~  (2) 

where  Hmc  is  the  incident  magnetic  field,  f  is  the  position  vector  from  the  origin,  k  is  the 
free  space  wavenumber,  and  the  time  dependence  has  been  suppressed.  The  incident 
plane  wave  gives  rise  to  electric  currents  directed  in  both  the  (f)  and  t  directions  on  the 
surface  S  as  well  as  scattered  far  fields  in  the  Or  and  (j)r  directions. 

3  MAGNETIC-FIELD  INTEGRAL  EQUATION 
3.1  Formulation 

Figure  3  illustrates  an  arbitrarily  shaped  perfectly  conducting  body  with  a  plane  wave 
incident  on  the  body.  By  definition  of  the  scattered  magnetic  field, 

H{f)  =  Hsc  +  Hin,.  (3) 

For  the  position  vector  f  not  on  the  surface  S 
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Figure  2.  Plane  Wave  Scattering  by  a  Conducting  Body  of  Revolution. 


Hsc{r)  =  -V  X  Asc{f) 

lio 

where 


Substituting  for  Asc{r)  in  Eq.  (5),  the  scattered  magnetic  field  becomes 
H.Ar)  =  ^  j  J(T')xV'*dS' 
where 

p-jk\T-f'\ 

■0  =  -p — =7r- 
\r  -  r  \ 

For  f  inside  the  sphere,  H{r)  equals  zero,  and  thus 
H,,{t)  =  -Hincir)- 

Substituting  Hsc{"f^)  from  Eq.  (8)  into  Eq.  (6)  we  obtain 


(4) 

(5) 

(6) 

(7) 

(8) 
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Figure  3.  Geometry  of  Perfect  Conductor  for  Plane  Wave  Scattering. 


-i?. 


(r )  =  ^  [  J{f')  X  V'rpdS'  (r  inside  S). 

47r  Js 

As  r  approaches  S  from  inside  the  surface,  Eq.  (9)  becomes^ 


(9) 


h  X  HUf)  =  ^  ^  X  |(f -  f')  x  /(r')l<iS'  (10) 

for  r  on  S,  where  we  have  used  the  vector  relationship 


j{f')  X  VV  =  -  f)  X  J{r').  (11) 

The  outward  normal  h  is  defined  as  h  =  u^y.ut.  Using  the  notation  in  Reference  6 
and  letting 

J{r')  =  <f>')  +  00  (12) 

where  u[  and  v!^  are  unit  vectors  in  the  t'  and  (j)'  directions  respectively,  Eq.  (10)  expands 
to 
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+  ^  / P'dt'  d4>'GJ\t\(l>'  +  mip'  -  p)cosv' 

—  (z'  —  z)  sint;']  cos(j)'  —  2/)C0ST;'sin^(Y)} 


+  —  [ p'dt'  [  d(f)'GJ^{t',  ((>'  +  <j))iz'  -  z)  sin  (j)' 
47r  J  Jo 

^  d<p'GJ\t',  (f)' +<!>){  y  sin  u  cos  u' 


p  sin  v'  cosv-{z'  -  z)  sin  v  sin  v']  sin  (f)'} 


/p£.Ti 

cl)'dt' J  #'GJ^t',0'+0){[(p'-p) 


(z'  —  z)  sin  u]  cos  (f)'  +  2p'  cos  v  sin^ ) } 


--  p)  cosv 


Tl  X  Hir] 


where 


_  1  +  jkR  jkR  (14) 


and 


In  Eq.  (13),  both  h  and  kt  are  to  be  evaluated  on  S  at  t  and  (f).  The  angle  between  Ut  and 
the  z-axis,  v,  is  positive  when  Ut  points  away  from  the  z-axis  and  is  negative  when  Ut  points 
toward  the  z-axis.  The  primed  variables  are  functions  of  t'  and  therefore  are  evaluated  at 

t'. 

We  want  to  take  advantage  of  the  geometry  of  the  body_^of  revolution  by  separating 
the  (f)  and  t  dependence  of  Jt  and  J<t,.  To  accomplish  this,  Jt  and  are  expanded  in  a 
Fourier  series  in  (j): 


OO 

n=— OO 

OO 


(16) 

(17) 


n=-oo 
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^' + ,f)  =  E 


The  Ut  and  u^,  components  of  Eq.  (13)  will  be  evaluated  using  the  orthogonal  properties 
of  sinusoids  and  the  following  integral  definition  of  Bessel  functions. 


^jx  COS  <i>~j 


for  m  =  0,  ±1,  ±2, . . . 


Rewriting  Eq.  (13),  we  have  for  the  Ut  component  of  the  incident  plane  wave 


2  4jr 


/r2Tr 

dt'p'fnit')  jf  d0V'”'^'G{[(p'-p)cosi;' 


—  {z^  —  z)  sin?;']  cos0'  —  2pcos?;'sin^(““)} 


J  dt'p'j’lit')  I  d(j)'e^’"^'G{z'  —  z)  sm4>' 


=  -- cos0te^*^^'°®^‘j"[J„4,i(A;psm0t)  +  Jn- lik  p  sin  St)] 


and  the  component 


m  , 

2  47r 


jdt'p'm 


d<j)'e^'^'^  G  [p'  sin  v  cos  v' 


p sin v' cos V  —  {z'  —  z)  sinv sin v']  sin (f)'] 


J  dt'p'jtit') 


d<f>'e^’''^'G{[(p'  -  p)  cos  V 


—  {z'  —  z)  sinu]  coB(j)'  +  2p'cost;sin^(— )} 

=  ^‘  {  27r  cos  v  sin  dtj  Jn  {kp  sin  $t)  — 

sinv cosdt‘Kj'^'''^[Jn+i{kpsin6t)  -  J„+i(fepsin^t)]}. 
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For  the  H‘^  incident  wave  we  have  for  the  iit  component 


^  +  ^  f  dt'p'jlit')  [  d<f)'e^^‘^'G{[{p' -  p)cosv' 

2  47r  J  Jo 

—  {z'  —  z)  sinu']  cos  4)'  —  2pcos'i;'sin^(^)} 

+  !f-  f  dt'p'jtit')  r  d(f>'e^'^^’G{z'  -  z)sin(t>' 

47r  J  Jo 

=  -  Jn-iikpsinOt)] 

2 

and  for  the  component 


2  47r 


—  psinu'cosu  —  (z' 


d(j)'e^'^‘^G[p'  sin  v  cos  v' 
—  z)  sin  V  sin  v']  sin  (j)' 


+  ^  j dt'p'jUt')  d4>'e‘''*'G{Hp'  -  p)cosv 
—  (2'  —  2)  sin  v]  cos  0'+  2p'cosusin^{^)} 


(23) 


(24) 


=  ^ sm»e'‘'^'”*'[;”'*-^.+i(*=P8inft)  -/J„-i(fc()sin«,)]. 

2 

The  t'  and  0'  integrations  can  be  converted  to  direct  summations.  With  the  t'  integral 
divided  into  discrete  segments,  the  value  of  the  integrand  is  evaluated  at  the  center  of  each 
segment  and  then  summed  over  the  number  of  f'-segments.  The  (f)'  integral  is  divided  into 
small  patches,  again  approximating  the  integrand  at  the  center  of  each  patch.  In  moment 
method  terminology,  this  ecjuates  to  using  pulse  basis  functions  and  unit  impulse  testing 
functions.  The  well-behaved  condition  of  the  magnetic-field  integral  equation  allows  such 

simple  basis  and  testing  functions  to  be  used. 

Rewriting  Eqs.  (21)  and  (22)  with  the  t'  and  (f)'  integrations  represented  as  summa- 

tions,  we  have 


M 


M 


^  +  At  ^  K“(t;,  tm)3im  +  ZJ  t^m  =  Ot) 

for  /  =  1,2,... 


m=l 

m^l 


m=l 


(26) 


and 


8 


(26) 


M  M 

+  At  ^  iff  (tl,  f  +  At  =  Sj(f,,  e.) 

m=l  m—l 

m^l  m^l 

for  1  =  1,2,..  .M 

where  n  is  the  mode  number,  I  is  the  number  of  the  t-segment,  and  m  the  segment  number 
for  the  t'  integration.  The  first  superscript  indicates  the  direction  of  the  current  and  the 
second  superscript  indicates  the  vector  component.  M  is  the  maximum  number  of  t-  and 
t'-segments,  and  At  is  the  length  of  the  t-segment. 

For  a  given  t{\)  and  t'{m)  the  only  unknown  quantities  in  Eqs.  (25)  and  (26)  are 
and  j^.  Thus  for  a  given  n,  t,  and  t',  Eqs.  (25)  and  (26)  can  be  represented  as 

-Ki\  ...  ifiv  if,'?  ...  Kit,  ]ri!i  rsr 

T^tt  T^tt  T^i<l>  -t  at 

^Ml  *  •  •  ^MM  *  ‘  •  ^MM  JM 

^11  ^11  •••  ^\M 

L^MI  •••  ^MM  ^M\  ^MM 

where  K  and  S  are  known  functions  of  t,  t',  (p',  and  t  and  6t,  respectively  and  the  /s  are 
the  unknown  Fourier  series  coefficients  of  the  current. 

The  superscripts  of  the  j  and  S  variables  indicate  the  vector  component  and  the 
subscripts  indicate  the  segment.  This  square  matrix  equation  can  be  solved  for  and 
for  each  n,  and  the  currents  can  be  found  from  y*  and  by  using  Eqs.  (16)  to  (19). 

Evaluating  the  <p'  integration  in  Eqs.  (21)  to  (24)  dX  t  =  t'  {p  =  p'  and  z  =  z')  yields 
R  =  2psin(072).  Therefore,  at  t  =  t',  the  phi  integration  does  not  converge.  This  lack  of 
convergence  is  not  due  to  the  magnetic-field  integral  equation  but  to  our  separation  of 
the  4>'  and  t'  integrations.  To  obtain  convergence  at  t  =  t',  one  must  evaluate  the  (p't'  area 
of  integration  more  carefully.  In  Reference  7,  it  is  shown  that  the  <p'  integration  can  be 
done  properly,  when  the  <j)'  integration  is  separated  from  the  t'  integration,  if  the  modified 
Green’s  function 


(28) 

replaces  G  in  Eq.  (14),  where 

R,=  yl(Atl2Y-Kp4/f 

(29) 

R,=,j(AtnYn^(pi/f 

(30) 

and  At  is  the  size  of  the  t-segment. 

3.2  Far  Field 

To  obtain  the  scattered  magnetic  far  field,  let  r  »  r'  in  Eq.  (6)  to  get 

Noting  that 

r  X  uf--cos{4’  -  +  cos#sin(i^  -  4'')ut 

r  X  U(  =  sin{^  -  0')  sin  +  cos  ^  cos{0  -  0')  sin  v'u^ 
fr-f'  =  f!  sin  e,  cos(0  -  0')  +  z'  cos  » 


(31) 

(32) 

(33) 

(34) 


allows  Eq.  (31)  to  be  written  as 


=  r  fi'dt'  r’d0'{j*(0'.«')i-cos(0-0')i 

47rr  Jo  Jo 

+  /(0'.(')|sin(0-  0')sln»1}e'‘l<’'""“<*-«+''“‘'' 


= 


-jke- 


jkr  pT 


pi  p^7V 

/  p'dt'  /  #'{  t')  [-  sinicj)  -  00  cos  9] 

Jo  Jo 


t')  [-  sin  0  cos  v'  +  cos(0  -  00  ^os  9  sin  v'] } 


^3k\p* s\nO  cos{<l}-4>')+ z' cos  6] 


where,  according  to  Eqs.,  (16)  and  (17) 


(35) 


(36) 


C50 

/(0'.*')  =  E 

n=-oo 


n=—oo 


(37) 

(38) 


The  superscripts  «  and  0  in  Eqs.  (35)  and  (36)  designate  the  conrponent  of  the  scattered 
field. 

Substituting  Eqs.  (35)  and  (36)  into  Eq.  (31)  and  using  Eq.  (20),  Eq.  (35)  is  rewritten 
as 
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(39) 


*'0  n=~oo 

[J„+i{kp'  sm0)  +  J„_i(A;p'sin^)]| 


and  Eq.  (36)  as 

(f)  =  r  p'dt' 

47rr  7o  n=-oo 
[J„+i(A:p'sin0)  -  J„_i(A:p'sin0)] 

sin  9  cos  v'Jn(kp'  sin  9)  —  cos  9  sin  v'j 

[Jn+i{kp'sm9)  -  J„_i(Vsin6>)]}). 

The  radar  cross  section  a  is  defined  as 

a  =  lim  47rr  - 

r^oo  n.  12 

■‘■■‘■tncl 


(40) 


(41) 


For  a  He  incident  wave  Eq.  (41)  becomes 


a  =  lim  47rr^ 

r-^oo 


(42) 


and  for  a  incident  wave  Eq.  (41)  becomes 


a  —  lim  47rr^ 

r— foo 


(43) 


where  the  first  and  second  superscripts  designate  the  components  of  the  incident  and 
scattered  fields,  respectively. 

For  the  plane  wave  incident  along  the  positive  z-axis,  only  the  ±  1  modes  are  needed 
in  the  Fourier  series  to  represent  J*  and  J'*’.  For  an  incident  plane  wave  off  the  z-axis,  the 
maximum  number  of  modes  Nmax  was  chosen  as 


Nmax  =  ( 1 . 1  fcr  sin  -f-  2)  (44) 

where  r  is  the  maximum  radius  of  the  body  of  revolution  measured  firom  the  z-axis.® 

®Yaghjian,  A.D.  (1977)  Near  Field  Antenna  Measurements  on  a  Cylindrical  Surface:  A  Source  Scattering  Matrix 
Formulation,  NBS  Technical  Note  696. 
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Figure  4.  Geometry  of  Perfect  Conductor  for  Dual-Surface  Magnetic-Field  Integral 

Equation. 


Typically  five  to  ten  segments  per  wavelength  are  used  for  both  the  number  of  t-segments 
and  phi-segments.  The  criterion  for  the  number  of  phi-segments  is  based  on  having  square 
patches  for  the  phi  integration.  Other  shaped  patches  were  used,  but  the  results  proved 
most  accurate  for  square  or  approximately  square  patches. 

4  DUAL-SURFACE  MAGNETIC-FIELD  INTEGRAL  EQUATION 
4.1  Formulation 

The  following  derivation  of  the  DMFIE  closely  follows  the  work  on  the  DMFIE 
presented  by  Woodworth  and  Yaghjian.^  The  DMFIE  can  be  derived  by  beginning  with 
Eq.  (10),  the  “interior”  or  “extended”  MFIE,  which  is  rewritten  here 

-Hinc{f)  =  ^  f  J{f')  xV'ipdS'  (f  inside  S).  (45) 

47r 

The  current  J(f)  in  Eq.  (45)  is  uniquely  determined  at  every  firequency  if  Eq.  (45)  is 
satisfied  for  all  r  inside  S.® 

Figure  4  shows  the  surface  S  and  a  surface  Ss  parallel  to,  and  some  distance  ^  >  0, 
inside  of  S. 

Waterman,  P.C.  (1965)  Matrix  formulation  of  electromagnetic  scattering,  Proc.  IEEE,  53:805-812. 
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By  adding  an  cross  Eq.  (64)  at  points  f  on  Ss  to  the  corresponding  points  on  S  in 
the  original  MFIE,  one  obtains  the  DMFIE 


n  X  Ho{f)  =  -n  X  ^  (f  J{f')  x  V  Vo(r,  f')  dS' 
2  47r  Js 


(46) 


where 


.ffo(fO  =  +  aHincir-  8n) 

(47) 

■00  (r,  f')  =  0(r,  f')  +  a:0(r  —  f') 

(48) 

as  usual 

1 

1 

'^1 

(49) 

II 

1 

The  DMFIE  will  have  a  unique  solution  for  J  at  all  real  frequencies  if  a  is  imaginary  and 
the  positive  real  constant  8  is  less  than  approximately  A/2.  For  the  body  of  revolution 
solution,  a  was  chosen  to  be  j  (\/^),  and  8  was  selected  to  equal  A/4. 

The  existing  MFIE  solution  is  easily  modified  to  obtain  the  DMFIE.  Hmc  and  xl)  are 
now  calculated  at  both  the  outer  surface  and  an  inner  surface  some  distance  8  inside.  This 
requires  that  computer  program  be  modified  to  simply  evaluate  these  functions  at  f—8h 
as  well  as  at  f. 

Figure  5  plots  the  normalized  radar  cross  section,  a/'Ka^,  in  the  backscattering 
direction  for  a  perfect  electric  conducting  sphere  as  computed  using  the  MFIE,  the  Mie 
series,  and  the  DMFIE.  The  spurious  resonances  in  the  MFIE  solution  are  eliminated  by 
the  DMFIE,  which  agrees  well  with  the  Mie  series. 


5  FFT  AND  CONJUGATE  GRADIENT  IMPLEMENTATION 
5.1  FFT  Method 

As  shown  in  Eqs.  (21)  to  (24)  the  expansion  of  the  current  into  a  Fourier  series 
allows  the  decoupling  of  the  (f)'  and  t'  integrals  in  the  integral  equation,  and  allows  the  <p' 
integration  to  be  performed  for  each  term  in  the  Fourier  series.  The  (j)'  integration  in  Eqs. 
(21)  to  (24)  can  be  rewritten  as 

p27r 

I{n)=  /  A(t,t',(/)')Ge^"^'#'  (50) 

Jo 

where  A{t,  t',  4>')  is  the  portion  of  the  equation  independent  of  n.  Eq.  (50)  has  the  form  of  a 
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Figure  5.  Backscattering  Cross  Section  Versus  ka  of  a  PEC  Sphere  Computed  with  the 

Mie  Series,  MFIE,  and  DMFIE. 


Fourier  transform  and  can  be  computed  efficiently  using  a  FFT.  lin)  must  be  computed 
at  a  discrete  integer  spacing  of  1,  (i,e.  n  =  1, 2, 3, . . .  Nj^sx)  according  to  the  criteria  in  Eq. 
(44),  and  so  the  FFT  must  be  sampled  at  discrete  intervals  A0'  =  27r/m  where  m  is  an 
integer.  This  allows  the  <f)'  integrations  for  all  of  the  modes  to  be  calculated  simultaneously. 
The  use  of  the  FFT  then  reduces  the  computational  time  required  to  perform  the 
integrations  but  increases  the  storage  required  to  store  the  additional  information.  If  the 
Fourier  series  requires  M.  modes  for  convergence  of  the  current,  it  is  necessary  to  solve  M 
systems  of  linear  equations  each  with  2N  unknowns. 

5.2  Conjugate  Gradient  Method 

The  conjugate  gradient  iterative  method  is  an  efficient  alternative  to  Gaussian  elim¬ 
ination  or  LU  decomposition  for  solving  the  large  systems  of  linear  equations,  AX  =  V- 
A  version  of  the  CG  method,  referred  to  as  Case  A  by  Sarkar  and  Arvas,^^  applies  to  a 
general  square  matrix  A  that  is  not  defined  to  be  symmetric,  positive  or  negative  definite. 
They  start  with  an  initial  guess  Xo  and  define 


7^o  =  V  -  AXo 


(51) 


Gedney,  S.D.  and  Mittra,  R.  (1990)  The  use  of  the  FFT  for  the  efficient  solution  of  the  problem  of  electromagnetic 
scattering  by  a  body  of  revolution,  IEEE  Ufans.  AntenTids  PTopagat.^  AP-38. 313-322.  \  r  u 

”  Sarkar,  T.K.,  and  Arvas,  E.  (1985)  On  a  class  of  finite  step  iterative  methods  (conjugate  dir^U^s)  for  the 
solution  of  an  operator  equation  arising  in  electromagnetics,  IEEE  Trans.  Antennas  Propagat.,  AP-22. 1058- 1066. 
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==  ^0  —  *4*7^0 


(52) 


and  iterate  for  i  =  1, 2, . . . : 


(53) 

Xi  =  Xi-\  +  aiPi-i 

(54) 

II 

1 

1 

! 

(55) 

Qi  =  A^Hi 

(56) 

(57) 

—  Gi  ^  biPi-i 

(58) 

where  ||^o||^  =  |cjP  is  the  inner  product  of  a  vector  with  itself  and  \cj\  denotes  the 
magnitude  of  the  complex  vector  element  Cj.  The  symbols  X,y^7l^Q,  and  V  are  vectors, 
a  and  h  are  scalars,  and  A*  denotes  the  Hermitian  conjugate  of  the  matrix  A.  The  CG 
method  minimizes  where  Hi  is  the  residual  vector,  and  gives  a  sequence  of  least 

square  solutions  (Xi)  to  AX  =  y  such  that  the  estimate  Xi  gets  closer  to  X  with  each 
iteration.  Iteration  continues  until  the  normalized  residual  error  (e  =  l|7?^t||/||!y||)  is  reduced 
to  a  chosen  value.  For  the  stated  problem,  X  is  the  vector  of  the  unknown  surface  currents, 
y  is  the  vector  of  the  incident  field,  and  A  matrix  represents  the  operator  on  the  right-hand 
side  of  Eq.  (46). 

The  sequence  Xi  will  converge  to  the  exact  solution  A"  in  a  finite  number  of  steps 
for  any  initial  guess  so  long  as  the  matrix  is  not  singular  and  the  round-off  errors  are 
negligible.  An  initial  guess  of  Aq  =  0,  the  null  vector,  was  used.  The  use  of  the  null  vector 
allows  the  flexibility  of  stopping  the  iteration  and  restarting,  using  the  current  estimate 
of  the  solution  as  the  new  initial  guess.  Restarting  eliminates  round-off  errors  but  requires 
an  additional  matrix  multiplication,  thus  slowing  convergence  and  increasing  the  solution 
run  time. 

The  radar  cross  section  of  two  right  circular  cylinders  with  equal  diameter  and  length 
was  computed  using  two  different  values  of  e,  10“^  and  10“®.  The  results,  in  Figures  6  and 
7,  show  no  difference  in  the  radar  cross  section  and  so  an  c  =  10“^  was  chosen  as  the  cutoff 
value  of  the  normalized  residual. 

As  stated  previously,  the  CG  method  in  practice  converges  more  rapidly  than  the 
theoretical  limit  of  iterations.  To  demonstrate  this  for  a  body  of  revolution,  the  total 
number  of  iterations  was  divided  by  the  total  number  of  modes  to  obtain  an  average 
number  of  iterations  for  the  problem  solution.  Figure  8  shows  the  slight  increase  in  the 
average  number  of  iterations  required  to  compute  the  bistatic  radar  cross  section  of  a  right 
circular  cylinder  with  the  same  diameter  and  length.  This  means  that  the  total  number 
of  iterations  required  increases  approximately  linearly  with  electrical  size  of  the  cylinder. 
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Figure  6.  Bistatic  Radar  Cross  Section  of  a  5A  by  5A  Right  Circular  Cylinder  with 
Normalized  Residual  Error  =  10  ^  and  10  ;  =  90  . 


The  CG  method  was  found  to  converge  steadily  for  the  well  conditioned  DMFIE  body  of 
revolution  solution  as  seen  in  Figure  9. 

5.3  Numerical  Implementation 

For  axial  incidence  of  the  plane  wave  {9t  =  0  or  180°),  only  the  ±1  Fourier  modes 
are  required.  For  an  incidence  angle  of  0°,  the  bistatic  scattering  of  a  300A  x  300A  right 
circular  cylinder  (6000  unknowns)  was  computed  using  Gaussian  elimination  to  solve  the 
set  of  linear  equations  for  five  segments  per  wavelength.  The  solution  required  570  minutes 
of  CPU  time  on  the  Cray  COO  supercomputer  and  is  shown  in  Figure  10.  The  size  of  the 
body  was  limited  by  the  COO  queue  restrictions,  not  the  available  RAM. 

Off  axis,  the  maximum  number  of  Fourier  modes  m  is  given  by  Eq.  (44).  The  number 
of  unknowns  N  is  proportional  to  d/ A,  the  electrical  arc  length  along  the  generating  curve. 
Using  a  right  circular  cylinder  of  equal  length  and  width,  for  each  m  the  matrix  fill  and 
solution  times  using  Gaussian  elimination  are  proportional  to  (d/X)?  For  broadside  inci¬ 
dence  m  is  proportional  to  d/X  so  the  total  matrix  and  fill  solution  times  are  proportional 
to  {d/Xy.  Using  the  FFT  implementation  discussed  previously  reduces  the  matrix  fill  time 
from  a  (d/A)'*  to  a  {d/Xflog2(d/X)  functional  dependence.  The  use  of  the  CG  iterative 
method  allows  one  to  reduce  the  (d/A)'*  solution  time  to  a  {d/Xf  solution  time.  This  means 
that  the  total  computer  time  required  for  the  body  of  revolution  solution  is  asymptotically, 
for  large  d/A,  proportional  to  {d/Xf  for  axial  incidence  and  (d/A)^  log2(d/A)  for  broadside 
incidence.  The  FFT  requires  an  increased  amount  of  computer  storage,  but  this  is  the 
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Figure  7.  Bistatic  Radar  Cross  Section  of  a  lOA  by  lOA  Right  Circular  Cylinder  with 
Normalized  Residual  Error  =  10"^  and  10“®;  0t  =  9O°. 


trade-off  one  must  tolerate  to  save  computer  time.  The  bistatic  scattering  of  a  20A  x  20A 
right  circular  cylinder  was  computed  at  broadside  incidence  {6t  =  90  degrees).  The  size  of 
the  cylinder  was  limited  by  the  available  RAM  on  the  Cray  COO.  The  solution  (Figure  11) 
required  10  minutes  of  CPU  time  for  Gaussian  elimination  with  no  FFT,  3  minutes  for 
Gaussian  elimination  with  the  FFT  and  approximately  2  minutes  for  the  CG  method 
with  the  FFT.  Figure  12  shows  the  CPU  time  versus  the  cylinder  radius  for  right  circular 
cylinders  with  the  same  diameter  and  length  using  the  FFT  to  decrease  fill-time,  CG  to 
decrease  solution  time,  and  using  strictly  Gaussian  elimination  to  reduce  matrix  solution 
time.  The  reduction  of  total  solution  time  of  electrically  large  bodies  by  the  FFT  and  CG 
method  is  clearly  demonstrated. 
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Figure  8.  Comparison  of  the  Average  Number  of  Iterations  Required  for  Solution  of  a 
Right  Circular  Cylinder  for  Five  and  Seven  Divisions  per  Wavelength;  Ot  =  90°. 


Figure  9.  Condition  Number  versus  Radius  of  a  Right  Circular  Cylinder  for  Five  and 
Seven  Divisions  per  Wavelength;  6t  =  90°. 
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Figure  10.  H-Plane  Bistatic  Radar  Cross  Section  of  a  300A  by  300A  Right  Circular 

Cylinder  with  e  =  10“^;  Ot  =  90°. 


BISTATIC  ANGLE  (deg) 

Figure  11.  H-Plane  Bistatic  Radar  Cross  Section  of  a  20A  by  20A  Right  Circular 
Cylinder  with  e  =  10“^;  0t  =  90°. 
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Figure  12.  Comparison  of  CPU  Times  for  FFT  and  CG  Method  and  Gaussian 
Elimination  for  Right  Circular  Cylinders;  0t  =  90°. 
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6  MODIFICATION  OF  THE  DMFIE  FOR  BODIES  WITH  SHARP  TIPS 


6.1  Formulation 

It  is  well  documented  that  the  original  MFIE  and  thus  the  DMFIE  degenerate  toward 
an  underdetermined  system  of  linear  equations  for  infinitesimal  thin  bodies.  5.12,13,14,15,16 
Reference  15  describes  the  numerical  implementation  of  the  MFIE  for  an  extremely  thin 
body,  a  metallic  disk  with  a  thickness  of  10~®A.  Reference  5  compares  the  measured  and 
MFIE-computed  values  of  the  E-  and  H-plane  RCS’s  for  a  perfectly  conducting  circular 
disk  with  a  radius  of  2A  and  a  thickness  of  A/20.  Reference  14  states  that  the  MFIE  is  not 
applicable  for  thin  or  open  bodies  but  does  not  define  what  the  criterion  is  for  a  thin  body. 
The  claim  that  the  MFIE  is  “useless  for  infinitely  thin  scatterers  but  finds  its  greatest  use 
in  large  smooth  conductors”  is  found  in  Reference  12. 

Consider  the  RCS  of  the  conesphere  shown  in  Figure  13  computed  using  the  DMFIE. 
Figures  14  and  15  show  E-plane  and  H-plane  monostatic  RCS  versus  target  aspect  angle 
for  the  DMFIE  and  measurements.^^  Note  the  large  discrepancy  between  the  E-plane 
computed  results  and  the  measurements  in  the  angular  region  from  approximately  0  —  60°, 
the  illumination  of  the  sphere  portion  of  the  conesphere.  After  measurement  error  was  ruled 
out  by  the  comparison  of  measurements  from  two  different  measurement  facilities,^®  the 
number  of  divisions  per  wavelength  in  the  f  and  (f)  coordinates  were  increased  to  ensure 
convergence  of  the  numerical  scheme.  Figure  16  compares  the  E-plane  DMFIE  numerical 
results  for  20,  30,  and  80  divisions  per  wavelength.  One  sees  the  extremely  slow  convergence 
of  numerical  solution  to  the  measured  results. 

The  E-  and  H-plane  DMFIE  RCS  results  for  the  conesphere  with  20  segments  per 
wavelength  were  compared  to  the  CFIE,  MFIE,  and  EFIE  results  from  CICERO,^®  a 
CFIE  formulation  for  axially  inhomogeneous  bodies  of  revolution.  The  CICERO  code  uses 
triangular  instead  of  pulse  functions,  which  may  account  for  the  slight  difference  in  the 
MFIE  and  DMFIE  E-pIane  solutions.  However,  Figure  18  clearly  demonstrates  that  the 
discrepancy  is  a  result  of  the  MFIE  formulation,  not  the  implementation  of  the  DMFIE. 
The  close  agreement  between  the  CFIE  results  and  measurements  suggests  that  the  CFIE 
is  the  “correct”  answer.  The  CFIE  solution  with  400  divisions  per  wavelength  will  be  used 
as  the  “exact”  reference  solution. 


'^Poggio,  A.J.,  and  Miller,  E.K.  (1987)  Solution  of  Three-Dimensional  Scattering  Problems,  in  Computer 
Techniques  for  Electromagnetics,  Mittra,  R.,  ed.,  Springer- Verlag,  p.  168. 

i^Sadasiv,  M.R.,  Wilton,  D.R.,  and  Glisson,  A.W.  (1982)  Electromagnetic  Scattering  by  Surfaces  of  Arbitrary 
Shape,  IEEE  Trans.  Antennas  Propagat.,  AP-3O:409-418. 

Glisson,  A.W. , and  Wilton,  D.R.,  (1980)  Simple  and  Efficient  Methods  for  Problems  of  Electromagnetic  Radiation 
and  Scattering  from  Surfaces,  IEEE  Trans.  Antennas  Propagat.,  AP-2 8: 593-603. 

Sancer,  M.I.,  McClary,  R.L.,  and  Glover,  K.J.  (1990)  Electromagnetic  Computation  using  Parametric  Geometry, 
Electromagnetics,  10:85-103. 

Phone  discussion  with  D.  Wilton,  University  of  Houston,  Houston,  TX. 

Larose,  C.  and  Shantnu,  M.,  of  David  Florida  Laboratory,  Canadian  Space  Agency,  Ottawa,  Ontario,  performed 
the  RCS  measurements.  One  measurement  was  performed  from  0  -  360°  in  aspect  angle,  but  the  symmetry  of  the 
conesphere  allows  the  measurement  data  from  180  -  360°  to  be  plotted  from  0  —  180°  aspect  angle 

Wing,  R.  and  Cote,  M.,  of  Rome  Laboratory,  Hanscom  AFB,  MA,  performed  a  second  set  of  measurements  on 
the  conesphere. 

Putnam,  J.M.  and  Medgyesi-Mitschang,  L.N.  (1987)  Combined  Field  Integral  Equation  Formulation  for  Axially 
Inhomogeneous  Bodies  of  Revolution,  McDonnell  Douglas  Corporation  Report  No.  QA003,  Volume  I:  Final  Technical 
Report,  St.  Louis,  MO.,  Sandia  National  Laboratory  Contract  No.  33-4257. 
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Figure  13.  Geometry  for  Conesphere. 


The  amplitude  and  phase  of  the  scattered  field  for  the  first  three  modes  were 
calculated,  individually  and  collectively,  over  our  aspect  angular  range  of  0  -  80°  for  the 
MFIE  and  CFIE  using  CICERO.  Figures  20  to  29  show  the  comparison  of  these  results 
for  the  different  formulations.^®  It  was  determined  firom  these  results  that  the  zero  order 
current  mode  is  not  properly  modeled  in  the  present  DMFIE  and  MFIE  formulations.  This 
assumption  was  further  verified  by  calculating  the  phase  and  amplitude  of  the  currents 
for  each  mode  at  an  incident  angle  of  20°  with  the  CFIE  solution.  Figure  30  shows  that  at 
this  incident  angle  only  the  zero  and  first  order  modes  are  significant  contributors  to  the 
total  current.  The  amplitude  and  phase  of  these  two  modes  are  shown  in  Figures  31-34  for 
the  DMFIE,  MFIE,  and  CFIE. 

Reference  15  states  that  the  cells  on  the  opposite  faces  of  a  thin  structure  can 
become  so  close  that  the  calculation  of  the  matrix  element  containing  the  cells  can  break 
down  numerically  and  give  erroneous  answers  to  the  MFIE.  The  presented  numerical 
solution  in  Reference  15  is  a  finer  integration  of  these  cells.  This  information  prompted 
the  modification  of  the  DMFIE  program  to  allow  finer  integration  of  the  self  ring,  that 
is  when  t  =  t',  for  the  zero  order  mode,  t-directed  current.  The  self  ring  patch  is  divided 
into  a  predetermined  number  of  tiny  segments  that  are  individually  integrated  and  then 
coherently  summed  to  obtain  the  self  ring  value.  Figure  35  shows  that  finer  integration  of 
the  self  ring  has  little  effect  on  the  DMFIE  solution.  The  only  visible  difference  in  the  five 
results  is  located  at  the  deep  null  at  approximately  15°. 


^“This  approach  was  suggested  by  J.M.  Putnam  of  McDoimell  Douglas  Research  Laboratories,  St.  Louis,  MO. 
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Figure  14.  Comparison  of  the  E-Plane  Monostatic  RCS  of  Conesphere  for  the  DMFIE 

and  Measurements;  ka  =  6.36. 


Singularities  of  the  tip  of  a  circular  cone  have  been  documented  by  various 
authors.^^’^^’^^  Reference  21  concentrates  on  the  circular  cone  shown  in  Figure  Al  and 
is  directly  applicable  to  the  conesphere  problem.  The  static  electromagnetic  expression, 
explained  in  additional  detail  in  Appendix  A,  was  implemented  for  the  self  ring  when  the 
source  point  was  at  the  tip  patch  and  the  field  point  up  to  twenty  patches  away  on  the 
cone  portion  of  the  conesphere.  Figures  36-38  show  the  current  and  RCS,  respectively,  at 
an  aspect  angle  of  20°  for  the  implementation  of  this  method  and  the  CFIE.  Figure  38 
shows  no  considerable  improvement  in  the  RCS. 

The  next  attempt  at  solving  the  problem  was  to  use  the  exact  tip  current  for  a 
semi-infinite  cone.^"*  From  Reference  24,  it  was  determined  that  the  t— directed  zero  order 
mode  current  at  the  tip  of  the  cone  has  a  Spherical  Bessel  function  ju{t)  variation  where 
t  is  the  radial  distance  from  the  tip  of  the  cone  in  wavelengths.  The  mathematical  and 
numerical  details  are  presented  in  Appendix  B.  The  t— directed  current  for  mode  zero  is 
now  approximated  by 


J“  =  Dj,{t) 


(59) 


Van  Bladel,  J.  (1983)  Field  singularities  at  the  tip  of  a  cone,  Proc.  IEEE,  71:901-902. 

^DeSmedt,  R,  and  Van  Bladel,  J.  (1986)  Fields  at  the  tip  of  a  metallic  cone  of  arbitrary  cross  section,  IEEE 
Trans.  Antennas  Propagat,  34:865-870. 

Jackson,  J.D.  (1975)  Classical  Electrodynamics,  Wiley,  New  York,  2nd  ed.,  pp.  94^102. 

^'*Trott,  K.D.  (1986)  A  High  Frequency  Analysis  of  Electromagnetic  Plane  Wave  Scattering  by  a  Fully  Illuminated 
Perfectly  Conducting  Semi-Infinite  Cone,  Dissertation  at  Ohio  State  University,  Columbus  ,Ohio. 
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Figure  15.  Comparison  of  the  H-Plane  Monostatic  RCS  of  Conesphere  for  the  DMFIE 

and  Measurements;  ka  =  6.36. 

where  D  is  some  constant.  The  program  was  then  modified  to  use  a  Spherical  Bessel 
function  variation  for  the  t-directed  current  for  the  segment  at  the  tip  of  the  conesphere 
with  a  segment  length  that  corresponds  to  the  location  of  the  first  maximum  of  the 
Spherical  Bessel  function.  Figures  39  and  40  show  the  comparison  of  the  current  fi:om  the 
modified  DMFIE,  MFIE,  and  CFIE.  One  sees  the  improvement  in  the  current.  Figure  41 
compares  the  original  DMFIE,  the  modified  DMFIE,  the  CFIE,  and  measurements.  One 
now  sees  that  the  modified  DMFIE,  CFIE  and  measurements  compare  fairly  well  in  the 
problem  region. 

We  find  that  the  failure  to  accurately  model  the  current  near  the  tip  affects  not 
only  the  current  near  the  tip,  but  also  many  wavelengths  away  fi-om  the  tip.  Therefore, 
the  expression  in  Eq.  (31)  for  the  scattered  magnetic  field,  which  shows  that  the  current 
within  a  quarter  wavelength  of  the  tip  has  little  effect  on  the  RCS,  does  not  account  for 
enough  of  the  physics  of  the  problem.  The  conclusion  is  that  for  the  conesphere  where  the 
t— directed  current  near  the  tip  for  mode  zero  has  a  strong  effect  for  many  wavelengths 
away  from  the  tip,  and  thus  a  strong  effect  on  the  RCS,  accurate  modeling  of  the  current 
must  be  performed. 
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Figure  16.  Comparison  of  the  E-Plane  DMFIE  Monostatic  RCS  of  Conesphere  for  20, 
50,  and  80  Divisions  per  Wavelength  with  Measurements;  ka  —  6.36. 
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Figure  17.  H-Plane  Monostatic  RCS  of  Conesphere  for  the  DMFIE,  MFIE,  EFIE,  CFIE, 

and  Measurements;  ka  =  6.36. 
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Figure  18.  E-Plane  Monostatic  RCS  of  Conesphere  for  the  DMFIE,  MFIE,  EFIE,  CFIE, 

and  Measurements;  ka  =  6.36. 


Figure  19.  E-Plane  Monostatic  RCS  of  Conesphere  for  the  MFIE,  CFIE,  and 
Measurements  with  400  t— Divisions  per  Wavelength;  ka  =  6.36. 
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Figure  20.  Amplitude  of  the  E-Plane  Scattered  Field  of  Conesphere  for  the  MFIE  and 

CFIE;  ka  =  6.36,  mode  zero. 


Figure  21.  Phase  of  the  E-Plane  Scattered  Field  of  Conesphere  for  the  MFIE  and  CFIE; 

ka  =  6.36,  mode  zero. 
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Figure  22  Amplitude  of  the  E-Plane  Scattered  Field  of  Conesphere  for  the  MFIE  and 

CFIE;  ka  =  6.36,  mode  one. 


igure  23.  Phase  of  the  E-Plane  Scattered  Field  of  Conesphere  for  the  MFIE  and  CFIE; 

ka  =  6.36,  mode  one. 


Figure  24.  Amplitude  of  the  E-Plane  Scattered  Field  of  Conesphere  for  the  MFIE  and 

CFIE;  ka  =  6.36,  mode  two. 


Figure  25.  Phase  of  the  E-Plane  Scattered  Field  of  Conesphere  for  the  MFIE  and  CFIE; 

ka  =  6.36,  mode  two. 
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Figure  26.  Amplitude  of  the  E-Plane  Scattered  Field  of  Conesphere  for  the  MFIE  and 

CFIE;  ka  =  6.36,  modes  zero  and  one. 


Figure  27.  Phase  of  the  E-Plane  Scattered  Field  of  Conesphere  for  the  MFIE  and  CFIE; 

ka  =  6.36,  modes  zero  and  one. 
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Figure  28.  Amplitude  of  the  E-Plane  Scattered  Field  of  Conesphere  for  the  MFIE  and 
CFIE;  ka  =  6.36,  modes  zero,  one,  and  two. 


Figure  29.  Phase  of  the  E-Plane  Scattered  Field  of  Conesphere  for  the  MFIE  and  CFIE; 

ka  =  6.36,  modes  zero,  one,  and  two. 
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Figure  30.  Amplitude  of  the  E-Plane  CFIE  Solution  Current  for  Conesphere;  ka  =  6.36, 

dt  =  20°,  all  modes. 


Figure  31.  Amplitude  of  the  E-Plane  DMFIE,  MFIE,  and  CFIE  Solution  Current  for 
Conesphere;  ka  =  6.36,  dt  =  20°,  mode  zero. 
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Figure  32.  Phase  of  the  E-Plane  DMFIE,  MFIE,  and  CFIE  Solution  Currents  for 
Conesphere;  ka  =  6.36,  $t  —  20°,  mode  zero. 
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Figvne  33.  Amplitude  of  the  E-Plane  DMFIE,  MFIE,  and  CFIE  Solution  Currents  for 
Conesphere;  ka  =  6.36,  9t  =  20°,  mode  one. 
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Figure  34.  Phase  of  the  E-Plane  DMFIE,  MFIE,  and  CFIE  Solution  Currents  for 
Conesphere;  ka  =  6.36,  9t  =  20°,  mode  one. 


Figure  35.  E-Plane  Monostatic  RCS  with  DMFIE  Using  Finer  Integration  on  Self-Ring 
for  t-Directed  Current,  Original  DMFIE,  and  CFIE;  ka  =  6.36. 
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Figure  36.  Amplitude  of  the  E-Plane  Modified  DMFIE  using  t''  Current  Dependence, 
Currents  of  MFIE  and  CFIE  Solutions  for  Conesphere;  ka  =  6.36,  6t  =  20°,  mode  zero. 


Figure  37.  Phase  of  the  E-Plane  Modified  DMFIE  using  t''  Current  Dependence, 
Currents  of  MFIE  and  CFIE  Solutions  for  Conesphere;  ka  =  6.36,  6t  =  20°,  mode  zero. 
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Figure  38.  E-Plane  Monostatic  RCS  for  Conesphere  for  the  Modified  DMFIE  using  t'' 
Current  Dependence,  Original  DMFIE,  and  CFIE;  ka  =  6.36. 


Figure  39.  Amplitude  of  the  E-Plane  Modified  DMFIE  using  Current  Dependence, 
Currents  for  MFIE  and  CFIE  Solutions  for  Conesphere;  ka  =  6.36,  dt  =  20°,  mode  zero. 
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Figure  40.  Phase  of  the  E-Plane  Modified  DMFIE  using  Current  Dependence, 
Currents  for  MFIE  and  CFIE  Solutions  for  Conesphere;  ka  =  6.36,  9t  =  20°,  mode  zero. 


Figure  41.  E-Plane  Monostatic  RCS  for  Conesphere  for  the  Modified  DMFIE  using 
Current  Dependence,  Original  DMFIE,  and  CFIE;  ka  =  6.36. 
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Figure  42.  Geometry  for  Right  Circular  Cone 


6.2  Numerical  Implementation 

The  original  DMFIE  code  was  modified  to  include  the  Spherical  Bessel  function 
dependence  of  the  tip  current.  The  code  was  automated  so  that  the  operator  can  at  will, 
use  this  Spherical  Bessel  function  calculation  of  the  better  tip  current  for  tipped  scatterers. 
A  subroutine  was  included  to  find  the  maximum  of  the  Spherical  Bessel  function  and  set 
the  length  of  the  first  segment  on  the  tip  of  the  scatterer  to  this  electrical  length.  This 
program  was  then  used  to  compute  the  RCS  of  conespheres  and  right  circular  cones  of 
various  values  of  ka. 

Figure  42  shows  the  geometry  and  dimensions  for  the  flat-backed  right  circular  cone. 
In  Figures  43-66,  the  RCS  versus  aspect  angle  for  different  electrical  sized  conespheres  as 
well  as  flat-backed  cones  is  presented.  The  CFIE  results  were  computed  with  CICERO 
using  100  divisions  per  wavelength  and  the  measurements  performed  as  stated  in  Reference 
17.  The  original  and  modified  DMFIE  results  used  20  divisions  per  wavelength. 

The  results  demonstrate  that  the  flat-back  cones  seem  to  be  less  sensitive  to  the 
current  at  the  tip  of  the  cone  than  the  conespheres  for  small  cone  half-angles.  There  is  no 
difference  between  the  two  DMFIE  results  when  the  half-cone  angle  is  greater  than  15°,  as 
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Figure  43.  E-Plane  Monostatic  RCS  for  Conesphere  for  the  Modified  DMFIE  using  ji, 
Current  Dependence,  Original  DMFIE,  and  Measurements;  ka  =  5.38. 


is  expected.  The  conespheres,  as  demonstrated  previously,  are  extremely  sensitive  to  the 
modeling  of  the  mode  zero  current. 
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ASPECT  ANGLE  (deg) 

Figure  44.  H-Plane  Monostatic  RCS  for  Conesphere  for  the  Modified  DMFIE  using  ju 
Current  Dependence,  Original  DMFIE,  and  Measurements;  ka  =  5.38. 


Figure  45.  E-Plane  Monostatic  RCS  for  Conesphere  for  the  Modified  DMFIE  using 
Current  Dependence,  Original  DMFIE,  and  Measurements;  ka  =  5.87. 
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Figure  46.  H-Plane  Monostatic  RCS  for  Conesphere  for  the  Modified  DMFIE  using  > 
Current  Dependence,  Original  DMFIE,  and  Measurements;  ka  =  b.87. 


ASPECT  ANGLE  (deg) 

Figure  47.  E-Plane  Monostatic  RCS  for  Conesphere  for  the  Modified  DMFIE  using  > 
Current  Dependence,  Original  DMFIE,  and  Measurements;  ka  —  6.86. 
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Figure  48.  H-Plane  Monostatic  RCS  for  Conesphere  for  the  Modified  DMFIE  using  > 
Current  Dependence,  Original  DMFIE,  and  Measurements;  ka  =  6.86. 


Figure  49  E-Plane  Monostatic  RCS  for  Flat-Backed  Cone  for  the  Modified  DMFIE 
using  Current  Dependence,  Original  DMFIE,  and  Measurements;  a  =  1.18A,  I  =  4.33A, 

a  =  15.26°. 
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Figure  50.  H-Plane  Monostatic  RCS  for  Flat-Backed  Cone  for  the  Modified  DMFIE 
using  >  Current  Dependence,  Original  DMFIE,  and  Measurements;  a  =  1.18A,  I  =  4.33A, 

a  =  15.26°. 


Figure  51.  E-Plane  Monostatic  RCS  for  Flat-Backed  Cone  for  the  Modified  DMFIE 
using  jv  Current  Dependence,  Original  DMFIE,  and  Measurements;  a  =  1.29A,  /  =  4.72A, 

a  - 15.26°. 
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Figure  52.  H-Plane  Monostatic  RCS  for  Flat-Backed  Cone  for  the  Modified  DMFIE 
using  >  Current  Dependence,  Original  DMFIE,  and  Measurements;  a=  1.29A,  /  =  4.72A, 

a  =  15.26°. 
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Figure  53.  E-Plane  Monostatic  RCS  for  Flat-Backed  Cone  for  the  Modified  DMFIE 
usine  Current  Dependence,  Original  DMFIE,  and  Measurements;  a=  1.39A,  I  =  5.11A, 

a  =  15.26°. 
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Figure  54.  H-Plane  Monostatic  RCS  for  Flat-Backed  Cone  for  the  Modified  DMFIE 
using  Ov  Current  Dependence,  Original  DMFIE,  and  Measurements;  a  =  1.39A,  I  =  5.11A, 

a  =  15.26°. 


Figure  55.  E-Plane  Monostatic  RCS  for  Flat-Backed  Cone  for  the  Modified  DMFIE 
using  u  Current  Dependence,  Original  DMFIE,  and  Measurements;  a  =  1.5A,  /  =  5.51  A, 

a  =  15.26°. 
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Figure  56.  H-Plane  Monostatic  RCS  for  Flat-Backed  Cone  for  the  Modified  DMFIE 
using  >  Current  Dependence,  Original  DMFIE,  and  Measurements;  a  =  1.5A,  I  =  5.51A, 

a  =  15.26°. 


Figure  57.  E-Plane  Monostatic  RCS  for  Flat-Backed  Cone  for  the  Modified  DMFIE  using 
Current  Dependence,  Original  DMFIE,  and  CFIE;  a  —  0.808A,  /  =  3.042A,  o:  15  . 
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Figure  58.  H-Plane  Monostatic  RCS  for  Flat-Backed  Cone  for  the  Modified  DMFIE  using 
ji,  Current  Dependence,  Original  DMFIE,  and  CFIE;  a  =  0.808A,  I  =  3.042A,  a  =  15°. 


ASPECT  ANGLE  (deg) 

Figure  59.  E-Plane  Monostatic  RCS  for  Flat-Backed  Cone  for  the  Modified  DMFIE  using 
Current  Dependence,  Original  DMFIE,  and  CFIE;  a  =  1.342A,  I  =  4.921A,  a  =  16°. 
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Figure  60.  H-Plane  Monostatic  RCS  for  Flat-Backed  Cone  for  the  Modified  DMFIE  using 
Current  Dependence,  Original  DMFIE,  and  CFIE;  a  =  1.342 A,  /==4.921A,  q:  =  16  . 


Figure  61.  E-Plane  Monostatic  RCS  for  Flat-Backed  Cone  for  the  Modified  DMFIE 
using  jiy  Current  Dependence,  Original  DMFIE,  and  CFIE;  a  =  1.71A,  I  =  5.0A,  a  =  20  . 
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Figure  62.  H-Plane  Monostatic  RCS  for  Flat-Backed  Cone  for  the  Modified  DMFIE 
using  jv  Current  Dependence,  Original  DMFIE,  and  CFIE;  a  =  1.71A,  I  =  5.0A,  a  =  20°. 
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Figure  63.  E-Plane  Monostatic  RCS  for  Flat-Backed  Cone  for  the  Modified  DMFIE 
using  jy  Current  Dependence,  Original  DMFIE,  and  CFIE;  a  =  2. 11  A,  1  =  5.0A,  a  =  25°. 
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Figure  64.  H-Plane  Monostatic  RCS  for  Flat-Backed  Cone  for  the  Modified  DMFIE 
using  ji/  Current  Dependence,  Original  DMFIE,  and  CFIE;  a  =  2.11A,  I  =  5.0A,  a  =  25  . 


Figure  65.  E-Plane  Monostatic  RCS  for  Flat-Backed  Cone  for  the  Modified  DMFIE 
using  ji,  Current  Dependence,  Original  DMFIE,  and  CFIE;  a  =  2.5A,  I  =  5.0A,  a  =  30°. 
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Figure  66.  H-Plane  Monostatic  RCS  for  Flat-Backed  Cone  for  the  Modified  DMFIE 
using  >  Current  Dependence,  Original  DMFIE,  and  CFIE;  a  =  2.5A,  I  =  5.0A,  a  =  30°. 
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7  CONCLUSIONS 


It  was  shown  that  the  DMFIE  eliminates  the  spurious  resonances  from  the  MFIE  for 
plane  wave  scattering  from  bodies  of  revolutions.  The  results  were  excellent  for  as  few  as 
seven  segments  per  wavelength  for  the  right  circular  cylinder,  using  simple  pulse  basis  and 


impulse  testing  functions.  w  u 

The  reduction  of  total  solution  time  of  electrically  large  bodies  of  revolution  by 

the  FFT  and  CG  method  is  clearly  demonstrated.  The  solution  times  using  Gaussian 
elimination  were  proportional  to  (d/X)^  for  axial  incidence  for  a  scatterer  with  equal  width 
and  length,  and  proportional  to  {d/Xf  for  broadside  incidence.  For  broadside  incidence,  the 
FFT  reduced  the  matrix  fill  time  from  a  {d/XY  to  a  {d/Xf  \og,{d/X)  functiona  dependence, 
and  the  use  of  the  CG  iterative  method  reduced  the  solution  time  from  (d/A)  to  (d/A)  . 
Therefore,  the  total  computer  time  required  for  electrically  large  bodies  of  revolution 
is  proportional  to  (d/Xf  for  axial  incidence  and  (d/Xflog^id/X)  for  broadside  incidence 
However,  use  of  the  FFT  increased  the  amount  of  computer  storage  required  from  a  [d/X) 


proportionality  to  a  (d/  A)^  proportionality. 

It  was  demonstrated  that  the  MFIE  and  the  original  DMFIE  do  not  correctly  produce 
the  dominant  current  for  scatterers  with  narrow  tips  and  therefore  gave  incorrect  results 
for  scatterers  with  narrow  tips.  By  inserting  the  exact  tip  current  near  the  tip  of  an  infinite 
cone  into  the  DMFIE  computer  program,  the  DMFIE  was  made  to  accurately  calculate 
the  scattering  from  bodies  with  narrow  tips.  This  allows  the  DMFIE  to  be  applied  to 
bodies  of  revolution  with  much  narrower  tips  than  the  original  DMFIE  solution  and  thus 
increases  the  number  of  scatterers  to  which  it  can  be  applied.  This  improved  solution  for 
cone  tips  has  been  incorporated  into  the  original  DMFIE  computer  program  m  a  way  that 


is  transparent  to  the  user.  „  , ,  .  ,  ,  •  r 

Our  future  plans  are  to  implement  the  dual-surface  electric-field  integral  equation  for 

perfectly  conducting  bodies -of  revolution.  The  dual-surface  electric-field  integral  equation 
would  serve  as  an  independent  check  on  the  solution  to  bodies  f ” 
allow  us  to  compute  the  RCS  from  infinitesimally  thin  conductors  for  which  the  DMFIE 


cannot  be  appUed. 
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Appendix  A 

Dominant  Term  of  Electrostatic  Tip  Current  Variation  for  the  Cone 


The  potential  is  expanded  in  spherical  harmonics, 


OO 

(f){R,  ^  AmuP^ {cos  d) 

Tn=0  V 


(Al) 


The  value  of  the  subscript  v  is  determined  from  the  requirement  that  4>  vanish  on  the 
cone,  which  leads  to  the  condition 


P™(cos^o)=0. 


(A2) 


Only  the  zero  order  mode  (m  =  0)  is  of  interest  as  this  is  the  mode  that  gives  the  incorrect 
current  value. 

Van  Bladel  in  Reference  21  presents  two  singularities,  the  electric  and  magnetic.  The 
electric  field  of  the  electric  singularity  is  not  bounded  as  is  the  magnetic  field.  The  outward 
normal  crossed  into  the  magnetic  field  results  in  a  current  in  the  f-direction  {ur  direction 
in  Reference  21).  The  radial  distance  from  the  tip  of  the  cone  is  t  {R  in  Reference  21). 
The  current  varies  as  t''  where  p  is  the  degree  of  the  Legendre  function  in  Eq.  ( A2) . 

The  right-hand  side  of  Eq.  (13)  for  the  Ut  component  of  J\  J“,  can  be  written  as 


(j)')  =  Jave 


(A3) 


where 

g={l+jkr), 


(A4) 
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"t"  ^Jave 

{^ty 


(A9) 

(AlO) 


Using  Eq.  (AlO)  in  Eq.  (A8)  one  obtains 

J“  =  / it'f  d4,’Bgt'r'^  (All) 

which  now  can  be  substituted  into  the  right-hand  side  of  Eq.  (13).  The  value  of  u  can  be 
determined  from  the  asymptotic  formula 


(A12) 
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Appendix  B 

Dominant  Term  of  Tip  Current  Variation  for  the  Semi-Infinite  Cone 


Starting  with  Eq.  (4.48)  in  Reference  24 


€mh‘'u\kd)  f  S{u,m,P)cosm(l>\j  ,  . 

^  [  p  I  T{u,  m,  P)  sin  m</)  J  ” 


Jp{kRo) 
^  kRo 


mP”  (cos  0o) 

(sin  6o) 


(Bl) 


^mJ-yi^Rp)  f  COS  m(f)  1 

^7(7-1- 5(7,  m,^)  sin  m(/)  j\  ^ 


-1- 


h^^\kd) 

kd 


] 


^P7(cos«„) 


where  the  geometry  for  the  semi- infinite  cone  and  the  application  of  a  point  source  are 
shown  in  Figures  Bl  and  B2.  The  radial  distance  from  the  cone  tip  is  Rq  and  the  incident 
angle  P  is  chosen  such  that  O<P<0o-=n-a,  where  a  is  the  cone  half-angle.  We  are 
solving  only  for  the  E-plane,  t— directed  current  for  mode  zero.  This  reduces  Eq.  (Bl)  to 


J!„W  = 


27r 


^  Q 


where 


T(7,0,/3) 


{P'‘{cos/S)(7  -  1)(7)  -  Pj(cos/3)} 

2 


(B2) 


(B3) 
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Figure  Bl.  Geometry  for  semi-infinite  cone.^^ 


ifw  — 


d^[P°(cos  0)]^sin^ 


g-7 


(B5) 


and  gamma  is  the  value  of  the  zeros  for  the  solution  of  J^(cos  Oo)  0-  The  source  will  be 
placed  in  the  far  field.  This  facilitates  the  use  of  the  large  argument  form  of  the  Hankel 
function  given  by 


lim  hl^^(x) 


(B6) 


and  therefore 

The  solution  is  now  specialized  to  the  case  of  plane  wave  incidence  by  moving  the  point 
source  far  away  and  letting 


(B8) 
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Figure  B2.  Geometry  for  application  of  dipole  moment.^^ 


Applying  this  to  Eq.  (B2) 

j$'  /Tj  \  _  ny  ( J7(fcfio)(27  +  l)f  {d/ d0\P^{cos pi)\)  \ 

7(7+ I)sin0o{d/d7[-f^(cos0o)]}  ) 


(B9) 


The  numerical  solution  of  Eq.  (B9)  requires  the  evaluation  of  Eq.  D.15  in  Reference  24,  a 
recursion  relationship  for  Legendre  polynomials,  at  m  =  0,6  —  P,  and  I'  —  j,  yielding 


d  m  (7  +  1)jP7+i(cos^)  -  (7  +  l)P7(cos^) 

- 


(BIO) 


Differentiating  Eq.(D.13)  in  Reference  24,  an  integral  equation  representation  for  the 
associated  Legendre  polynomial,  with  respect  to  gamma  and  evaluating  at  m  =  0,9  =  do, 
and  1/  =  7,  is  also  necessary  for  the  evaluation  of  Eq.  (B9)  and  resiilts  in 


^Pj{cos9o)  =  — 
07  TT 


-  sin[(7  +  .5)0;] 
(cosx  —  COsOY'^^ 


(Bll) 


The  substitution  of  Eqs.  (BIO)  and  (Bll)  into  Eq.  (B9)  facilitates  its  numerical  evaluation. 

The  first  zero  can  be  determined  fi:om  Eq.  (A12),  tables^®  or  computer  programs.  The 
third  approach  was  used  for  a  cone  half-angle  a  of  12.5°  resulting  in  uq  =  0.226,  ui  =  1.326, 

“Hall,  R.N.  (1949)  The  application  of  non- integral  Legendre  functions  to  potential  problems,  Journal  of  Applied 
Physics,  20:925-931. 
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Figure  B3.  Amplitude  of  the  E-Plane  Current  From  the  Semi-Infinite  Cone  Formulation 
for  Three  Zeros  and  CFIE  Solution  for  Conesphere  for  an  Arc  Length  of  Three 
Wavelengths;  ka  =  6.36,  9t  =  20°. 


and  1/2  =  2.415.  Figs.  B3  and  B4  show  the  current  that  results  using  the  first  three  values  of 
u  in  Eq.  (B9)  compared  with  the  CFIE  solution  of  the  conesphere  at  a  =  12.5  and  6t  =  20°. 
The  semi-infinite  solution  using  three  zeros  is  the  best  fit  to  the  CFIE  at  a  distance  of 
approximately  one-half  wavelength  from  the  tip  of  the  body.  It  would  seem  that  one  should 
use  this  solution,  but  further  investigation  shows  that  Eq.  (B9)  for  multiple  values  of  u  is 
a  function  of  cone  and  incident  angles  with  a  complex  phase.  This  requires  the  use  of  only 
one  u  in  solving  Eq.  (B9)  and  Figure  B4  shows  the  agreement  to  be  acceptable  to  the  first 
maximum  of  the  current. 


62 


0  0.05  0.1  0.15  0.2  0.25  0.3  0.35  0.4  0.45  0.5 
ARC  LENGTH  (wavelengths) 


Figure  B4.  Amplitude  of  the  E-Plane  Current  Prom  the  Semi-Infinite  Cone  Formulation 
for  Three  Zeros  and  CFIE  Solution  for  Conesphere  for  an  Arc  Length  of  One-Half 

Wavelength;  ka  =  6.36,  9t  =  20°. 
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Appendix  C 

First  Maximum  of  Spherical  Bessel  Function 


The  need  for  the  maximum  of  the  Spherical  Bessel  function  j^{z)  was  demonstrated 
in  Appendix  B.  There  are  a  number  of  available  schemes  to  calculate  the  value  of  this 
function.  We  chose  one  that  was  computationally  fast  and  easy  to  implement.  We  decided 
to  use  the  series  representation  of  where 


Mz)  (9)  ^  juir/,,  I  fc . 


^2'  f^kir{u  +  k  +  i) 


(Cl) 


and 


Jn{z)  = 

Taking  the  first  two  terms  of  the  series,  we  have 


(C2) 


j1 _ 1. 

u  +  2.5)  J 


(C3) 


After  some  algebraic  reduction 

r  ^  ^(*'+2)  "i 

dn{z)  =  v^l  2('^+i)r(i/ +  1.5)  ~  r(i/  +  2.5)2(‘'+3)  j  ■ 

The  derivative  of  jn{2)  is  defined  as 


(C4) 


Abramowitz,  M.  and  Stegun,  LA.  (1972)  Handbook  of  Mathematical  Functions,  NBS 


(C5) 


Jni^J  -  V  TT  j  2{>^+i)r(i^  +  1.5)  r(iy  +  2.5)2(‘^+3)  j 
Set  y„(z)  equal  to  zero  and  solve  for  z  to  get 


where  z  is  the  radial  distance  in  wavelengths  from  the  cone  tip  for  our  case. 
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